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Abstract 

We consider the equation ut = u^x + b{x)u{l — u), x M, where 
b{x) is a nonnegative measure on R that is periodic in x. In the case 
where b{x) is a smooth periodic function, it is known that there exists 
a travelhng wave with speed c for any c > c*(6), where c*{b) is a 
certain positive number depending on b. Such a travelhng wave is often 
called a "pulsating travelling wave" or a "periodic travelling wave", 
and c*{b) is called the "minimal speed". In this paper, we first extend 
this theory by showing the existence of the minimal speed c*{b) for 
any nonnegative measure b with period L. Next we study the question 
of maximizing c*{b) under the constraint J^^ j^s^b{x)dx = aL, where 
a is an arbitrarily given constant. This question is closely related 
to the problem studied by mathematical ecologists in late 1980's but 
its answer has not been known. We answer this question by proving 
that the maximum is attained by periodically arrayed Dirac's delta 
functions aL X^fcGZ ^{^ ^ kL). 
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1 Introduction 



Travelling wave solutions describe a wide class of phenomena in combus- 
tion physics, chemical kinetics, biology and other natural sciences. From the 
physical point of view, travelling waves usually describe transition processes. 
Transition from one equilibrium to another is a typical case, although more 
complicated situations may arise. Since the classical paper by Kolmogorov, 
Petrovsky and Piskunov in 1937, travelling wave solutions have been inten- 
sively studied. For example, the monograph of Volpert, Volpert and Volpert 
[12j provides a comprehensive discussion on this subject. 

From the ecological point of view, travelling waves typically describe the 
expansion of the territory of a certain species, including, in particular, the 
invasion of alien species in given habitat. This kind of process may occur in 
both homogeneous and heterogeneous media. Models for biological invasions 
in spatially periodic environments were first introduced by Shigesada et al. 
in dimensions 1 and 2 (see [9l [lOl E]). More precisely, they considered 
spatially segmented habitats where favorable and less favorable (or even 
unfavorable) zones appear alternately and analyzed how the pattern and 
scale of spatial fragmentation affect the speed of invasions. In their study, 
the spatial fragmentation was typically represented by step functions which 
take two different values periodically. Mathematically, their analysis was 
partly unrigorous as it relied on formal asymptotics of the travelling wave 
far away from the front. 

Berestycki, Hamel [1] and Berestycki, Hamel, Roques [5] extended and 
mathematically deepened the work of Shigesada et al. significantly, by deal- 
ing with much more general equations of the form = V • {{A{x)Vu)) + 
f{x,u) in with rather general smooth periodic coefficients and by devel- 
oping various mathematical techniques to study the effect of environmental 
fragmentation rigourously. 

Among other things, they proved that, under certain assumptions on 
the coefficients, there exists c* > such that the equation has a pulsating 
travelling wave solution if and only if c > c*. Furthermore, they showed that 
the minimal speed c* is characterized by the following formula: 

c* = min{c > | 3A > such that ^(c. A) = 0}, 

where //(c, A) is the principal eigenvalue of a certain elliptic operator associ- 
ated with the linearization of the travelling wave far away from the front. A 
more detailed account of this result will be stated in Subsection 12.11 in the 
special context of our problem. 

By using a totally different approach Weinberger [Tl] also proved the 
existence of the minimal speed c* of pulsating travelling waves in a more 
abstract framework. His method relies on the theory of monotone operators 
and is a generalization of his earlier work [13] to spatially periodic media. 
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It is important to note that, as far as one-dimensional diffusion equations 
are concerned, the minimal speed c* coincides with the so-called spreading 
speed for a large class of monostable nonlinearities. Here the "spreading 
speed" roughly means the asymptotic speed of an expanding front that starts 
from a compactly supported initial data (see Definition 12.81 for details). An 
early study of spreading speeds in multi-dimensional spaces can be found 
in [21 El [13]. Weinberger [14j then studied the spreading speeds of order- 
preserving monostable mappings and applied the results to spatially periodic 
reaction diffusion equations and lattice systems. 

Berestycki, Hamel and Nadirashvili [6j also studied the spreading speed 
of reaction diffusion equation in a very general periodically fragmented envi- 
ronment where both the coefficients and the domain itself are periodic. Re- 
cently, the same authors |7j studied the spreading speed of reaction-diffusion 
equations with constant coefficients, but in very general domains which are 
not necessarily periodic. 

In this paper we consider the following equation 

ut = Uxx + b{x)u{l — u) , X G M, (1) 

where b{x) is either a smooth function or a measure satisfying b{x) > and 
b{x + L) = b{x), X G R, for some L > 0. 

By the above-mentioned work [H O [6l [U], the minimal speed c* of 
travelling waves is well-defined at least as far as 6(x) is a smooth function, 
and it coincides with the spreading speed. We denote this minimal speed by 
c*{b). The goal of the present paper is to consider the variational problem 

Maximize c*(b) 
b 

under the constraint 

/ b{x)dx = aL, (2) 

J[0,L) 

where a > is an arbitrarily given constant. In other words we want to 
find out whether or not there exists an optimal b{x) that gives the fastest 
spreading speed. We will show that the maximum of c*{b) does indeed exist 
but that it is not attained by any smooth function b{x) but by a measure 
which is composed of periodically arrayed Dirac's delta functions. 

In order to study the above problem, we have to consider the equation 

ut = Uxx + b{x)u{l — u), X G M, (3) 

where 6 is a nonnegative measure satisfying (j35|) . We introdcue two impor- 
tant quantities c*(b) and c*(5). The former denotes the minimal speed of 
travelling waves for equation The latter is a quantity associated with 
a generalized eigenvalue problem (see Definition 12.131 below) . It has been 
known that c*{b) = c*(6) if b is smooth. As we will see later, c*{b) = c*(6) 
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even if 6 is a measure (Theorem 12.161 ) . We will then show in Theorem 12.171 
that the maximum of c* (b) is attained by 

b{x) = h{x) ■.= aL^6(^x- {k + ^)Ly (4) 

where 5{x) is Dirac's delta function. 
Change begins 

Strictly speaking, we have to distinguish the travelling wave speeds in 
the positive direction and those in the negative direction (see Definition 12. 6|) . 
The above mentioned quantities c*{b) and c*(6) are associated with travelling 
waves in the positive direction. However, as we will see later in Theorem 
12.161 the two speeds - positive and negative - are always equal, therefore no 
am biguity occurs by not specifying the direction of the travelling wave. 

Change ends 

This paper is organized as follows. In Section 2, we introduce basic no- 
tations and state the main results. In Section 3, we prove the well-posedness 
of equation ([3]). In Section 4, we consider a generalized eigenvalue problem 
associated with equation We show that c*(6) is bounded when b is 
smooth and that c*(6) is bounded when b varies in a certain set of measures. 
In Section 5, we show that the minimal speed c*{b) of travelling waves exists 
when 6 is a measure, and it coincides with c*(6) and also with the spread- 
ing speed. Then we complete the proof of Theorem 12.171 In Section 6, we 
prove the lemmas on equicontinuity of the solutions of Cauchy problem (see 
equation (fT^ below). These lemmas are used in Sections 3 and 5. 

2 Notation and main results 



2.1 Basic notation 

In this subsection we introduce some notation and recall some known results 
which will be used later. 

In what follows we fix constants L > and a > 0. Let A(a) be the set 
defined by 

A(a) := {6(x) G C^(M) | b{x) > 0, 6(x) = 6(x + L) and / b{x)dx = aL}. 

J[0,L) 

Definition 2.1. A(a) is defined to be the sequential closure of A(a) in the 
space of distribution on M. More precisely, b G A(a) if and only if there 
exists a sequence in A(a) such that 

/ b{x)ri{x)dx = lim / bn{x)ri{x)dx (5) 

for any test function rj G C^(M), where the left-hand side of (0) is a formal 
integration representing the dual product <b,T]>. In this sense, we denote 
that bn^b in the weak * sense. 
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Since each bn is positive, 6 is a nonnegative distribution. Consequently, 
5 is a Borel measure on M. Therefore, ^ holds for every 77 G Co(M). 

In what fohows, we wih not distinguish the measure b and its density 
function b{x), as long as there is no fear of confusion. Thus we wih often 
use expression as in the left-hand side of ([5]). 

We will also note that, since 6yj(x) is L — pcviodic^ is also L — pGviodic 
in the following sense: 

/ b{x)r]{x + L)dx = / b{x)ri{x)dx (6) 

for 77 G Co(M). 

The following lemmas will be useful later: 

Lemma 2.2. Let C A(a) he a sequence converging to some b G A(a) 
in the weak * sense. Let r/(x) he a continuous function on M satisfying 

00 

Emax \r](x + kL)\ < 00. (7) 

fe=— 00 

Then r] is b-integrahle on M and the following hold: 

lim / bn{x)ri{x)dx = i b{x)r]{x)dx, (8) 



„ 00 

/ b{x)\ri{x)\dx < aL > max \ri{x + kL)\. (9) 
Jr 0<x<L 
fe=— 00 

Proof. For each integer M > 0, we define a cut-off function quix) by 

{0 for l^l > (M + 1)L 

1 for |x| < ML 

M + 1 — L~^|x| otherwise. 

Then, since qM{x)\r]{x)\ is continous and compatctly supported, we have 

b{x)qM{x)\ri{x)\dx = lim / bn{x)qM{x)\ri{x)\dx. (10) 
: Jr 

Note also that, since each 6„ belongs to A(q), 

00 „i 

fen(a^)gM(a;)|r/(x)|(ix = / 6„(x)g'M(2; + A;L)|r/(x + A:L)|(ix 



A;=— 00 ' 
00 



< qL > max |r7(x + A:L)|. 



0<a;<L 
fe=— 00 
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This and 1^ imply 

» oo 

/ b{x)qM{x)\r]{x)\dx < aL > max \'q{x + kL)\ 
Jr , 0<x<L 



k=—oo 



for M = 1, 2, 3, • • • . Letting M — t- oo and applying the monotone conver- 
gence theorem, we obtain Q. Hence ry is 6-integrable on M. 
Next we observe that 

bn{x)r]{x)dx - / hn{x)qM{x)r]{x)dx < / bn{x){l - qM{x))\'n{x)\dx 
JR JR 

< aL y max Inix + kL)\ 
" ^ 0<x<L ' ^ ^' 

\k\>M ~ ~ 

The assertion ([8]) easily follows from this and (jlOp . The lemma is proven. □ 

Lemma 2.3. Let {bn} C A(a) be a sequence converging to some b £ A{a) 
in the weak * sense. Let rj{x) be an L-periodic continuous function on M. 
Then 



n— >-oo 



Jim / bn{x)r]{x)dx = / b{x)r]{x)dx, (11) 

'[0,L) J[0,L) 



Proof. By the bounded convergence theorem, we have 

b{x)ri{x)dx = lim / b{x)pe{x)r]{x)dx, (12) 



'[o,L) 

where Pe is a cut-off function defined by 



Pe{x) = < 



for X G (— oo, — e) U [L, oo) 
e~'^{x + e) for e,0) 

1 for X G [0, L - e) 
e~^{L — x) for xG[L — e,L). 



On the other hand, by the L-periodicity of 6„ and ??, we have 

bn{x)ri{x)dx = / bn{x)p£{x)ri{x)dx (n = 1, 2, 3, • • • ) 



'[0,L) 

for any < e < L/2. It follows that 



lim 

n—>oo 



[0,L) 



bn{x)r]{x)dx = lim / bn{x)pe{x)r]{x)dx = / b{x)pe{x)r]{x)dx. 



Combining this and (jl2p . we obtain the desired identity. The lemma is 
proven. □ 
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We remark that, because of the L— periodicity of the measure b £ A(a), 
we can regard b as functional in the space C*(]R/LZ), which we denote by 
[6], in the fohowing sense: 

[b]{x)[r]]{x)dx = / h{x)r]{x)dx 
R/LZ -'[0,L) 

for any [q] G C(M/-LZ). Here r]{x) is the periodic continuous function on M 
associated with \rj\. Furthermore, by (jlip . we have if 6„ — ^ ^ in the weak * 
sense, then — )• [6] in the weak * sense in C*(M/LZ). 
We consider the following reaction-diffusion equation: 

i^t = Uxx + b{x)u{l — u), X S M, 

where b G A(a). Here we have not specified the range of t, but what we have 
typically in mind is either t > or t G M. Obviously A(a) C A.(a), therefore 
([3]) is a generalization of (1). 

Definition 2.4. Let / C M 6e any open interval. A continuous function 
u{x, t) : M X / — )■ M is called a weak solution of JH]) for t £ I (or a solution 
in the weak sense) if for any rj{x,t) G C^(M x /), 



ur]t dxdt = I [urjxx + b{x)u{l — u)ri) dxdt, 
xi JrxI 

where the second integral on the right-hand side is understood in the following 
sense: 

/ b{x)u{l - u)r] dxdt = / (V / b{x)u{l - u)ri dx) dt . (13) 

JrxI J I ^f^J[kL,{k+l)L) ' 

We next consider the following Cauchy problem: 

{ut = Uxx + b{x)u(\ — u) (x G M, t > 0), 

u{x,0) = uo{x) >0 (xGM), 

where uq G C(M) n L° 



(14) 



' oo / 



Definition 2.5. A continuous function u{x, t) : M x (0, oo) — t- M is called a 
mild solution of (fH|) if 

limu{x,t) = Uo{x) for any x G M 
and if it can be written as 

u{x,t)= G{x - y,t)uo{y)dy 



+ / / G{x - y,t - s)b{y)u{y,s){l - u{y,s)) dyds, 
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where ^ 

G(x, t) := , exp ( — ^) . 

As we will show in Section 3, a mild solution of (jl4p always exists for 
any 5 S A{a) and G C(M) n L°°(M) with uq > 0, and it is unique and is 
a weak solution. 

It is easily seen that, if u{x,t) is a mild solution of (|14p . then for any 
constant r > 0, u{x,x + t) is a mild solution of (fM|l with initial data u{x, r) 
(see Remark l3.5p . 

We call a function u{x, on M x M a mild solution for t G M if, for any 
r £ M, u{x,t + r) is mild solution of ()14p with initial data no(x) = u{x,t). 

Definition 2.6. A mild solution u{x, t) of Q for t € M is called a trav- 
elling wave solution (in the positive direction) if < u{x,t) < 1 and if 
there exists a constant T > such that 

u{x- L,t) = u{x,t + T) /or (x,t) G M X M, 
lim u{x,t) = 1, lim u{x,t) = locally uniformly in t G M. 
/i is called a travelling wave solution (in the negative direction) if 

u{x + L,t) = u{x,t + T) /or (x,t) G M X M, 
lim n(x, t) = 0, lim n(x, t) = 1 locally uniformly in t G M. 



Change begins 



Remark 2.7. /n what follows, unless otherwise specified, by a travelling 
wave we usually mean the one in the positive direction. 



Change ends Here we call the quantity c := L/T the speed (or the 
average speed or the effective speed) of the travelling wave solution u{x,t). 

Berestycki and Hamel Berestycki, Hamel and Rogues [5] and Wein- 
berger established the existence of the minimal speed of travelling wave 
solutions for general monostable nonlinearities f{x, u) satisfying certain con- 
ditions, and they also gave an eigenvalue characterization of the minimal 
speed (see (fTSj) below). Here, for f{x,u) = b{x)u{l — u), if 6 G A(a), it is 
not difficult to see that f{x,u) satisfies the assumptions in [H [5l [H]. So 
for equation ([1]), we know that for any b G A(a), there exists the minimal 
travelling wave speed c*{b) > in the following sense: 

c > c* (6) =^ There exists travelling wave with speed c ; 
< c < c*(6) =^ No travelling wave with speed c exists. 

To be more precise, c*{b) is defined to be the minimal travelling wave speed 
in the positive direction. As mentioned in Remark 12.71 one can also define 
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the minimal travelling wave speed in the negative direction, which one may 
call c* (6) . As we will explain in Theorem 12.161 we always have c* (6) = c* (b) , 
even when b{x) is not symmetric, therefore we do not need to distinguish 
the two minimal wave speeds. 

As we have mentioned earlier, c*{b) also coincides with the so-called 
"spreading speed" of expanding fronts for ([3]). Here, we define the spreading 
speed as follows: 

Definition 2.8. A quantity c**{b) > is called the spreading speed (in the 
positive direction) if for any nonnegative initial data no ^ with compact 
support, the mild solution u{x, t, uq) of ([3]) satisfies that 

(i) lim u{x,t,uo) = uniformly in {x > ct} if c > c**(b), 

(ii) lim u{x,t,Uo) = 1 uniformly in {0 < x < ct} if < c < c**{b). 

A quantity c**{b) > is called the spreading speed (in the negative di- 
rection) if for any nonnegative initial data uq ^ with compact support, the 
mild solution u{x,t,UQ) of ([3]) satisfies that 

(i) lim u{x,t,uo) = uniformly in {x < —ct} if c> c**{b), 

(ii) lim u{x,t,Uo) = 1 uniformly in {—ct < x < 0} if < c < c**{b). 
Change begins 

Remark 2.9. Basic properties of the spreading speed in general periodic 
environments are studied in ]^\^\^ IM^ - It is known, at least for smooth b, 
that the spreading speed c**(b) coincides with the minimal wave speed c*{b) 
(and, similarly, c**{b) = c*(b)). As we will show later in Theorem \2.1(A we 
have c**{b) = c**{b) for any b G A(a). 

It is known that in the "leading edge", namely the area where u ~ 0, we 
have the asymptotic expression 

~ e-^(^-'=*)V(a;), (15) 

where tp{x + L) = ip{x) > 0, and A > is some constant. Substituting (fT5]) 
into equation ([T|), we obtain the identity 

- <(x) + 2XtP'{x) - b{x)ip{x) = (A^ - Ac)V'(x). (16) 

This observation motivates us to introduce the following operator, which 
generalizes the operator on the left-hand side of (|16p to the case where b{x) 
is a measure: 
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Definition 2.10. Forb £ A(a), we define an (unbounded) operator —L^^i on 
the Banach space {V' G C(R) | il){x) = ip{x + L)} with \\ip\\ = max^gK |V'(a^)| 
as follows: 

-L^li/jix) = -il^"{x) + 2Xil)'{x) - b{x)il){x). 

Here the derivatives are understood in the "weak sense" by which we mean 
that —L)^ iil) = g if and only if, for any ip G C^(M), 

/ (—93" — 2A93' — 699) ip dx = I ipg dx. 
Jr Jr 



Change ends 



Definition 2.11. For b G A(a), A > 0, w;e call b) the principal eigen- 
value of the operator —L-^^ i := —tp"(x) + 2Xip'{x) — b{x)^p{x) if there exists 
a positive continuous function ip with ip{x) = 'ip{x + L) satisfying 

- i;"{x) + 2XtP'{x) - b{x)^p{x) = fi{X, (17) 

in the weak sense. Here ip is called the principal eigenfunction. 

For the principal eigenvalue, the following proposition holds. We will 
prove it in Section [5] Change begins by converting the problem (jl7p into a 



more regular eigenvalue problem for a compact positive operator. Change ends 



Proposition 2.12. For any b G A{a) and A > 0, the principal eigenvalue 
/x(A, b) exists, and it is unique and simple. 



Change begins We also note that the principal eigenfunction of (|17p 



belongs to H^^^iW), as we will see in Subsection 4.2. Now observe that 
A^ - Ac in (fT6l) is a constant and that ip > 0. Therefore, ()16p implies that 
n{X, 6) = A^ — Ac if b is smooth. In view of this, we define the following 



quantities when 5 is a general measure. Change ends 



Definition 2.13. We define the minimal speed in the positive direction c*(b) 
and a related value c*(6) as follows: 

c*{b) := inf{c > | traveling wave in the positive direction with speed c exists} 
c*(6) := inf{c > | 3A > such that ^jl{X, b) = X^ - Ac}. 

Similarly, the minimal speed in the negative direction c*(b) and c*(6) as 

c*{b) := inf{c > | traveling wave in the negative direction with speed c exists} 
c*(6) := inf{c > | 3A > such that fi{-X, b) = X^ - Ac}. 
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Change begins 

If b{x) is a smooth nonnegative function, then by the results of [H [5l |6] 
or by those of jl4] . the foUowing proposition holds. 

Proposition 2.14 ([HEIE], [H]). For smooth b G A(q), we have 

c**ib) = c*ib) = clib), c**{b) = c*{b) = cm. (18) 

Change ends 

We will see later that conclusion of above proposition remain to hold 
when 6 is a measure. 

2.2 Main results 

We are now ready to present our main results. 

It is shown in [5] that for any b G A(a), we have c*(6) > 2y/a = c*(q). 
Our first theorem gives an upper bound on c*{b) : 

Proposition 2.15. For any b € A(a), 

2^/^<c*{b) < 2y^a + a^L^. (19) 
Moreover, ifb^a, then < c*{b). 

As we mentioned above, the inequality 2^/a < c*{b) is found in [5]. The 
main novelty of this theorem is the upper bound. 

While the lower bound in p9|) is sharp since the equality holds for b = a, 
it has not been known whether c*{b) attains its maximum in A(a) or not. 
The next two theorems shows that c*{b) does not attain its maximum in 
A(a) but it does in the extended class A(a). 

Theorem 2.16 (Minimal speed). For any measure b £ A(a), there exists 
c*(b) > such that a travelling wave in the positive direction with speed c 
exists if and only if c> c* (6) . In other words, c* (b) is the minimal travelling 
wave speed in the positive direction. Furthermore, 

c*{b) = c:{b). 

Similarly, the minimal travelling wave speed in the negative direction c*{b) 
exists and c*{b) = c*(6). Furthermore, for any b G A(a), 

c*{b) = c*{b). 

Theorem 2.17 (Optimal coefficient). 

c*(/i) = sup c*(6) = nrax c*(6). 

beA{a) beA(o) 

Moreover, 

c*{h) > c*{b) for any b G A(a). 
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Change begins 



Theorem 2.18 (Spreading speed). For b € A(a), the spreading speed in the 
positive direction c**(b) and that in the negative direction c**{b) exist and 

c**{b) = c*(b), c**(b) = c*(b). 

Consequently, 

c**{b) = c**(b). 



Change ends 



To prove Theorem 12. 171 the following proposition is important. We will 
prove it in Subsection 14.31 



Proposition 2.19. Let {b„} be a sequence in A(a) converging to some 
b G A(a) in the weak * sense. Then 

c*(6) = lim c*(6„). 



3 Reaction-diffusion equation with a Borel-measure 
coefficient 

In this section, we establish the well-posedness of the Cauchy problem (|14|) . 

Theorem 3.1. For any given nonnegative initial data uq £ C(M) n L°°(M), 
the problem p4|) has a unique mild solution. This mild solution is also 
a weak solution and it depends continuously on the initial data in the L°° 
norm. 

Proof. First we show that the solution u{x, t) exists in the weak sense. Let 
bn G A(a) satisfy 6„ — )• 5 in the weak * sense. Then for any given initial data 
uq(x), the problem 

f iUn)t = iUn)xx + 6„(x)u.„(l - U„) -^g) 
\ Un{x,0) = Uo{x) 

has a classical solution Un{x,t) for any n G N. By the comparison principle, 
< Un{x,t) < 1}. Hence u„ (n G N) are uniformly bounded. 

Let < ti < ^2 be two positive numbers. By Lemma 16.31 which we will 
prove in Section 6, the family of solutions {un{x,t),x G M, i G [ti,t2]} are 
uniformly equicontinuous with respect to x and t. Here the modulus of 
equicontinuity may depend on ti and ^2- Applying Arzela-Ascoli theorem, 
we can get a subsequence, which we still denote by {un{x, t)} that converges 
uniformly in {x,t) G [-M,M] x [ti,t2] for every M > and < ti < t2- 
The limit function u{x, t) = lim„_j.oo Un{x,t) is defined for every {x,t) G 
M X (0, +oo) and satisfies 

ut = Uxx + b{^)u{l - u) (21) 
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in the weak sense. Next we show that u satisfies 



limti(x,t) = uq{x) for any x G 



(22) 



To see this we first note that, by Lemma 13.21 below, u can be written in the 
form 

u{x,t)= G{x - y,t)uo{y)dy 

JR 

+ / / G{x - y,t - s)b{y)u{y,s){l - u{y,s))dyds. 
Jo 

For any x G R, the first integral 

/ G{x — y,t)uo{y)dy ^ uo{x), as t — )• 0. 
Jr 

By Lemma 12.21 the second integral can be estimated as follows 
G{x -y,t- s)b{y)u{y, s)(l - u{y, s)) dyds 



< 



D 



zds, 



where D is a constant depending on ||u(-, t)||/^oom\. Since ||n(-,t)|| roof'in>\ IS 
bounded, we have 



G{x — y,t — s)b{y)u{y, s)(l — u{y, s)) dyds — )■ 0, as t — )• 0. 



Consequently ([22]) holds. 

Next, let uo,uo G C{M.) n L°°(M) be arbitrary and let u,u be the cor- 
responding weak solutions of ()14p . the existence of which has been proven 
above. Then w := u — u satisfies 



Wt = Wxx + rn{x, t)w (x G M, t > 0), 
-u;(x,0) = wq{x) {x G M), 



(23) 



where m := b{x){l — u — u) \s a measure-valued function of t. By Lemma 
3.21 below, we can express w as 

w{x,t)= / G{x -y,t)wo{y)dy 

rt . (24) 
+ / / G{x-y,t-s)m{y,s)w{y,s)dyds. 



Define 



pit) = \\wi-,t)\\j 
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Then, since 

G(x — y,t — s)m(y, s)w(y, s) dy < \\m(-, s)\\ max G{x — y,t — s)\w(y, s)\, 
we have 

p{t)<p{0) + M f P^'}^' (25) 
Jo V47r(t - s) 

where M is a constant such that ||m(-,t)|| < M for t > 0. By Lemma 7.7 of 
Alfaro, Hilhorst, Matano [1], it follows that 



Consequently 

||n(-,t)-{i(-,t)|Uoo(K) <e^^*/^(l + -= / j^ds)\\uQ-uo\ 



This proves the continuous dependence on the initial data and the unique- 
ness of the mild solution. 

□ 

Lemma 3.2. The function u{x, t) in the proof of Theorem \3.1\ is also a mild 
solution. In other words, u{x,t) can be written as 



u{x,t) = G{x -y,t)uo{y)dy 

^ .t . _ (27) 

+ / / G{x-y,t- s)h{y)u{y,s){l-u{y,s))dyds, 
Jo Jr 

where the second integral on the right-hand side is understood as in (|13p . 

Proof. Denote u{y, s){l—u{y, s)) by f{u{y, s)). Since u„ is a classical solution 
of equation ()20p . it is also a mild solution. Hence 



Unix,t) = / G{x -y,t)uo{y) + / / G{x - y,t - s)hn{y)f{un{y, s)) dyds. 

JO 

To prove this lemma, it is sufficient to prove 







G{x -y,t- s)b{y)f{u{y, s)) dyds 







lim / / G{x - y,t - s)bn{y)fiuniy,s)) dyds. 



We consider 

ft 

G{x -y,t- s) [bn{y)f{un{y, s)) - b{y)f{u{y, s))) dyds, 
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which we denote by Zn{x,t). First we have 



Znix,t)\< / / G{x - y,t - s){bniy) -b{y))f{u{y,s))dydi 



G{x -y,t- s)bn{y) {f{u{y, s)) - f{un{y, s))) dydt 



(28) 



Since u{x, t) is bounded and {6„ — 6} are uniformly bounded hnear func- 
tionals on any interval [/cL, [k + 1)L] where G Z, by Lemma 12.21 we have 



G{x -y,t- s){bn{y) - b{y))f{u{y,s))dy 



< V 2aLmax G{x-y-kL,t-s)\\f{u)\\L^(^y,iot]) 



< 2aL||/(u)||ioo(iRx[o,t]) 
for some constant Ci > 0. Since 



Gi 



lo Vt- s 
by Lebesgue convergence theorem, 



ds < +00, 



lim 



G(x -y,t-s) {bn{y) - h{y)) {f{u{y, s))) dyds = 0. 



Similarly, 



Co 



G{x - y,t- s)bn{y){f{u{y,s)) - f {un{y , s))) dy < -j= 
for some constant G2 > 0. Again, by Lebesgue convergence theorem. 



lim 



G{x -y,t- s)bn{y) {f{u{y, s)) - f{un{y, s))) dyds = 0. 



Hence 



u{x,t)= G{x - y,t)uo{y)dy 

JR 



+ J I G{x -y,t- s)b{y)u{y,s){l-u{y,s))dyds. 



□ 
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Prom Lemma [6. II in Section 6 and Theorem l3.H the following proposition 
can be obtained easily. 

Proposition 3.3. Let u{x,t,Uo) be the mild solution of p4|) with ini- 
tial data uq. Then for any to > and M > 0, the family of functions 
^o)}||uo||lo°{k)<m uniformly equicontinuous in x. 

Consider the construction and the uniqueness of mild solution. The 
following comparison principle holds. 

Proposition 3.4 (Comparison principle). Let u{x,t,UQ) be as in Proposi- 
tion l^.^i Then uq < vq implies 

u{x, t, Uq) < u{x, t, Vq) for X G M, t > 0. 

Change begins 

Remark 3.5. If u{x,t) is a mild solution of (|14p . then for any constant 
r > 0, u{x,x-\-t) is a mild solution of ()14p with initial data u{x,t). This is 
obvious if b is a smooth function, since a mild solution is a classical solution. 
In the general case where b G A(a). We can take a sequence of smooth bn 
with bn ^b (in the sense of ??) and use the approximation argument found 
in Lemma \3.S\ 

Change ends 



4 The linear eigenvalue problem 
4.1 Basic estimates 

We recall that /^(A, b) denotes the principal eigenvalue of the problem 

-V'" + 2AV'' -b^ = fi{X, b)ip, ipix + L) = Tpix) 

(see Definition 12. lip . In this subsection, we estimate fi{X,b) both from the 
above and below. Here we introduce the following notation. Let 

/fpV := W e HlM) I ^{^) = + L)} 

with the norm 

J[0,L) 

Let 

= {V e H^er\i^{x) > 0}. 

Note that the following embedding is compact: 
H^^^ ^ C7(M)nL°°(M). 
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Lemma 4.1. For any b G A(q;) and any ip{x) G El, we have 

f 2p'^{x)dx- [ b{x)i)'^{x)dx>-{a + a^L^)f tp'^{x)dx. 

J[0,L) J[0,L) J[0,L) 

Proof. First, for any xi,X2 G [0,L\, we have 

if^{x2) - V'^(xi) = / 2'4){x)'4)'{x)dx. 
Hence, for any positive number > 0, we have 



if^{x2) - i^'^ixi) <\ I i)''^{x)dx + k I il?{x)dx. 

k .l[Q,L) J[0,L) 



Multiplying the above inequality by 6(^2) and integrating it by {xi,X2) G 
[0,L] X [0,L], we get 



L / b{x2)tp'^ {x2)dx2 — ceL / ip'^{xi)dxi 

J[Q,L) J[0,L) 

<aL^(\ I ij'^{x)dx + k I i;^{x)dx). 
^kJ^o,L) J\0,L) ) 



'[0,L) J%L) 

This is equivalent to 

L / b{x)ip'^ {x)dx — aL / ip'^{x)dx 
J[o,L) J[o,L) 

<aL^(\ I ij'^{x)dx + k [ i;^{x)dx). 

^kJ[0,L) J[0,L) ^ 

Letting A; = aL, we obtain 

/ tp'^{x)dx- [ b{x)'il)^{x)dx>-{a + a^L^) I ^l?{x)dx. 

J[0,L) J{Q,L) J[0,L) 

Lemma 4.2. For A > 0, 6 G A(a), it holds that /Lt(A, 6) > /x(0, 6). 
Proof. Let ^ > be the principal eigenfunctions of — LA,b- Then 

+ 2AV'' - = /u(A, 6)V'. 

Multiplying this by ip and integrating it from to L, we have 

/ ip''^{x)dx— / b{x)Tp'^{x)dx = ii{\,b) / ip'^{x)dx, 

J[0,L) J[0,L) J[0,L) 



□ 
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hence 



MA, 6) 



I[o,L) '^'^{x)dx - /[o^^) h{x)i^^{x)dx 



On the other hand, for ;u(0, 6), we have the variational formula 

. /[0,L) '^'^ (^)^^ - /[0,L) ^(^)'^^ (^)^^ 
/i(0,6) = mm — ^ ' ' . 

Consequently /x(A, 6) > fj.{0,b). 

Lemma 4.3. For any b G A(a), holds that 

— a > /i(A, 6) > —a — o?L^. 
Proof. Dividing (|17p by ■(/'(a^) and integrating it from to L, we get 

Ui\M+a)L+ [ (—)^dx = 0. 
J[o,L) ^ 

It implies that fi < —a. By Lemma l4.ll and 14.21 

/i(A, b) > /i(0, b)>-a- a^L^. 



(29) 
□ 

(30) 
(31) 



□ 



Lemma 4.4. There exists a constant F > such that for any A > 0, 
b G A(a), the principal eigenfunction ip S El of the operator —Lx^b satisfies 

max ih „ 
mmip 

Proof. By equation ([3T]l . for any xi,X2 G [0, L), 
^'^ — ' In V'(x2) — ln^(xi)| 



In 



il^{xi) 



X2 ^/ 



< 



/ |^|<ix<VZj/ {^fdx 

J[0,L) W V •^[O.i) ^ 



< LV-(M(A,6) + a). 
Then for any xi, X2 G [0, L], 



This implies 



Note that 



-LV-(MA,fe)+a) < V'(^l) < gV-(MA,fe)+a)_ 
~ il{x2) ~ 



^/-{n{X,b) + a) < VoPL^ = aL. 
Therefore, by setting F = e"^^ , we obtain the desired estimate. 



□ 
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4.2 Uniform bounds of c*(6) for smooth b 

In this subsection, we will prove Proposition 12.151 Before doing that, we 
recall the following pointwise max-min formula for the principal eigenvalue. 

Proposition 4.5. For b G A(a), 

-r{^) + 2X^'{x)-b{x)i;{x) 
uiA,b) = max mi — . 

i>eELnC'^{R) xm ip[x) 

We omit the proof of the above proposition as it is easy. A more gen- 
eral version of the above proposition is fonud in [ 4 ; Proposition 5.7]. The 
following lemma easily follows form the definition c* (6) in Definition 12.131 

Proposition 4.6. For b G A(a), 

. -l^{\,b) + X'^ 

c„(0) = mm ; . (32) 

A>o A 

Proof of Proposition \2.15[ By Proposition I2.14^ it suffices to show that 



2\/a < c*(6) < 2\/a + a2L2. 
This follows immediately from (|30p and ()32p . The proposition is proven. □ 



4.3 Uniform bounds of c*(6) when 6 is a measure 

In the previous subsection, we discussed the boundedness of c*{b) for b E 
A(a). In this subsection, we derive the same bounds for c*(6) when b S A(q;). 
We start with the fo llowing proposition: 



Change begins 



Proposition 4.7. Let bn be a sequence in A(a) converging to some b in the 
weak * sense and let A„ G M 6e a sequence converging to some A G M. Then 
there exists a constant 13 such that 

n{\n-, bn) — >• /5 as n ^ oo. 

Furthermore, 

P = fi{X,b). 

Proof. Lemma HTHl shows that {/u(A„, bn)} is uniformly bounded. There exist 
functions ^pn G El with HV'nllL"" = 1 satisfying 

-^p'^ + 2A„Vn - bntpn = /"(A, bn)tpn- 

By Lemma [4. 41 there exists a constant F > such that max'0„/min'0„ < F. 



Change begins Now we multiply 

+ 2An< - (/i(A, bn) + bn{x))iJr, 
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with ^„ and integrate it from to L. We get 



[0,L) J[0,L) 



Hence 



(V'^)^dx+ / {i;nfdx= / {l + ^{X,bn) +bn{x))'4'ldx. 

[0,L) J[0,L) J[0,L) 

By the fact that HV'nllL"" = 1 and (j30p . we see from the equation above 
that {llV'nll/fi^^} is uniformly bounded. Change ends Then there exists 
^p G Hper such that for some /3 there exists a subsequence {V^n^}, 

V'nfc ~^ ^ weakly in Hp^j. and strongly in C(M) 

and ^(Anj.,^nfc) — >• /3 as — )■ +00. Therefore 

—ip" + 2AV'' — bip = /3ijj in the weak sense, 

So /i(A, 6) = /3. We will prove the uniqueness of fj,{X, b) in Proposition 
15.121 Check Therefore, it is not difficult to see fj,{Xn,bn) — >• fj.{X,b). 

□ 



Change ends 



The following corollary follows immediately from Proposition 14.71 

Corollary 4.8. Let bn he a sequence in A(a) converging tob £ A(a) in the 
weak * sense. Then fi(X, bn) — )• ^(A, b) locally uniformly in X > 0. 

We also remark that the continuity of A 1— n{X, b) follows easily from 
Proposition 14.71 

Proof of Proposition \2.19l By Definition 12.131 we have 

c*(5„) := inf{c > 0| 3A > such that fJ.{X,bn) = X^ - Ac} 
c*(5) := inf{c > | 3A > such that ^(A, b) = X^ - Ac}. 

Since //(A, 6„) — )■ /i(A, b) locally uniformly in A and since /i(A, 6^) is uniformly 
bounded by Lemma 14.31 We immediately obtain the conclusion. 

□ 

Combining Proposition 12.151 and 12.191 a-iid recalling that every b £ A(a) 
can be expressed as a weak * limit of sequence in A(a) (see ([5])), we obtain 
the following proposition. 

Proposition 4.9. For any b G A(a;), 

2Va < cl(b) < 2\/~aVo?I?. 
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4.4 Maximizing c*(fe) 

In Subsection 14.31 we have shown that c*(6) is bounded. In this subsection, 
we consider the variational problem 

Maximize c* (b) 

6eA(a) 

and show that the maximum is attained by h{x) defined in (jj]). 
Lemma 4.10. 

c*(/i) = rnax c*(6). 

beA{a) 

To prove the above result, we need some auxiliary lemmas. 
Lemma 4.11. // there exist A > 0, c > and ip G El such that 

-< + 2AV''-(A2-Ac + /i(x))V' = 0, (33) 
then Tp is a piecewise function and ip^L/l) > ip{x) for x G M. 
Proof. First, we just integrate (p3]l from to L to get 

(A^ - Xc + h{x))^pdx = 0. 

[0,L) 

This implies A^ — Ac < 0. 

Since h{x) = for a; G (-L/2,L/2), 

-V" + 2AV'' - (A2 - Ac)^ = {-^<x<^) 

in the classical sense. Furthermore, ip is continuous up to x = ib-| since El C 
Consequently by the classical maximum principle and the positivity 
of ^{x), along with the negativity of A^ — Ac, the maximum of ip must be 
attained at x = ±^ . Hence 

i^{^)>Hx), for xG (-^^,^^). (34) 
From the periodicity of il^{x), we get that 

ipi^) > ipix), X G M. 
Since can be expressed as 

V;(x) = cie^i^ + cae'^^^' m(-^,^) 

for some constants ci, C2, z^i, f2, ~ 0) ^^'^ ''P'i^ + 0) both exist. Conse- 
quently, is piecewise on M. □ 
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Lemma 4.12. For any b G A(a), c*(/i) > c*{b). 

Proof. In what follows c will stand for c*(/i) for simplicity. Then, by the 
definition of cl{h), there exists ^ G El such that 

-Tp" + 2XiP' - (A^ - Ac + h{x))tlj = 

in the weak sense for some A > 0. We define the convolution b * ip for 
L— periodic functions by 



[b * 'iIj]{x) := iIj{x -y)b{y)dy. 

J[0,L) 

Then one can easily see that tp := b*i/j belongs to El n C^(M) since b is 
and ip is piecewise C^. Consequently, 

= 6 * ( - V" + 2A^' - (A^ - Ac + b{x))ip) 
= -b* ip" + 2Xb*Tp' -b* ((A^ - Ac + h{x))ip) 
= -ip" + 2X1/;' - (A^ - Ac + b{x))ip - (6 * {hiP) - b{b * ip)) 



-i>" + 2Xip' - (A^ - Ac + b{x))ip - {aLiP{^)b - b{b * ^)) . 



Note that 



aL'ilj{^)b{x) - bix)[b * ip]ix) =b{x){aLiP{^) -[b* ip]ix)) 

=b{x){[b*iP{^)]{x)-[b*iP]{x)) 

=b{x)[b*{i,{^)-^{x))]{x) 
>0 

by ()34p . Furthermore, since the left-hand side of the above inequality is 
continuous and periodic in x, we have aLip{^)b{x) — b{x)b * ip > Xa for 
some cr > 0. It follows that 

-Ip" + 2A'0' - (A^ - A(c -a)+ b[x))'p > 0, 

hence, by Proposition 14.51 

Combining this and Proposition 14. 6t we obtain 

c*(6) <c-a<c:= cl{h). 
Moreover, by Proposition 12.14] 

c*{b)<c:{h) 

for any b G A(a). □ 
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Proof of Lemma \4-10[ The proof is immediate from Lemma [4.12l and Propo- 
sition [2l9l □ 



5 Pulsating travelling waves 

In this section we wih prove our main results: Theorems 12.161 12.171 12.181 

In Subsection 5.1, we show that the minimal speeds c*(6) and c*(6) exist 
and that they coincide with the spreading speeds c**{b) and c**{b), respec- 
tively. In Subsection 5.2, we prove c*(6) = c*(5) and c*(5) = c*(6) thereby 
completing the proof of Theorem 12.161 In Subsection 5.3, we complete the 
proof of Theorems 12. 171 and 12. 181 Without loss of generality, in this section, 
let L = 1. 

5.1 The existence of the minimal speed and spreading speed 

Our strategy here is to use the general results of Weinberger [14J and Liang 
and Zhao [5]. For this purpose, we need to consider the solution semiflow 
of ([3]) on the space C(M; [0, 1]) with respect to the local uniform topology, 
where the convergence Un ^ u means that Un{x) — >• u{x) uniformly on 
any bounded interval. First, we consider the linear space i?C(M) = {u\u £ 
C(M)nL°°(R)} which contains C(M; [0, 1]). Equip BC{R) with the topology 
of locally uniform convergence. It is easy to show that this linear topological 
space i?C(M) has the following properties. 

Proposition 5.1. For each m £ Z and n G N, define 

max \u{z)\ 




2W 



Then each \\ ■ is a norm on BC{M) and for any M > 0, \\ ■ ||^ defines a 
topology equivalent to the local uniform topology on C(M; [—M,M]). 
Furthermore, for any u G i?C(IR), 
Wrong begins 



\\u\\jn < ll'^llm ^ "'ll'^llmi G ^, f^- G N 

and 

2"l"l||u||^ < ll^li;; < 2l™l||'u||;}, Vm e Z,n e N. 



Wrong ends 






Modification begins 



n 

u\\ln+j < \\u\\m < X] W'^W^+i G Z,n E N, j = 0, . . . ,n - 1 
1=1 
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and 

2"l™l||n||[J < < 2l"^l||'u||[J, Vm e Z,n G N. 



Modification ends 



Proposition 5.2. C(M; [0,1]) is a bounded dosed subset of BC{R). More- 
over, C(M; [0,1]) is complete. 



In Section 3, we have shown that the mild solutions of ()14p depend on 
the initial data continuously in L°°(R). In the following proposition, we will 
show that the continuous dependence also holds with respect to the local 
uniform topology. 

Proposition 5.3. The mild solutions of (|14p depend on the initial data 
continuously in the space C(M; [0, 1]) . Precisely, for any e > and T > 0, 
there is some rj > such that for any two solutions u, v of |j^[ ) with initial 
data uq, vq, if \\uo — vo\\l < rj, then i) — i)||o < e for t G [0, T]. 

Proof. Let w = u — v with wq = uq — vq, then 

w{x,t)= / G{x - y,t)wo{y)dy 

ft r _ (35) 

+ / / G{y,t — s)b{x — y)m{x — y, s)w{x — y, s) dy ds, 



Jo JR 

where m(x,s) = l—u{x,s)—v{x,s). By comparison principle, < uo{x),vo{x) < 
1, Vx S M implies that < u{x,t),v{x,t) < l,Vx G M, t > 0, and then there 
\m{x, s)| < 1. Let 

h{t,x) = / G{x -y,t)wo{y)dy 

and 

h{t,x)= / / G{y,t-s)b{x-y){l-u{x-y,s)-v{x-y,s))w{x-y,s)dy ds. 
Jo Jr 

First, we consider Ii. For any e > and T > there is some rj > 0, such 
that if \\wo\\l = ||/i(0, Olio < V, then ||/i(t, Olio < ^ any t G [0,r]. Then 
we consider I2. By Lemma 12.21 



G{y, t — s)h{x — y)m{x — y, s)w{x — y, s) dy 

<MaLy max G{y,t — s) max \w{x — y,s)\ 
^yeli,i+i] j/e[j,j+i] 

where M is a constant such that ||m(-, t)|| < M for t > 0.. 
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Hence, we have 



max |/2(t,x)| < MaL / max > max G{y,t — s) max \w{x — y,s)\ds 

x&[k,k+l\ Jq xe[k,k+l] ^ ye[i,i+l] ye[i,i+l] 

< MaL / > max G{y,t — s) max \w{z,s)\ds. 
Jo f^ye[i,i+i] ze[k-i-i,k~i+i] 

This implies that 

Whit, -Mil < MaL [ V max t - s) V 2"l^l max |it;(z, s)|)<is 

= MaL / V max Giy,t - s)\\wi-,s)f_i_^ds 

<MaL f y max t - s)2l*+2| ||^(.^ 

Jo f-^3/e[j,j+i] 



< 



t 



lo yjt- 

for some positive constant C . Then we have 

r-t Ql 

lk(-,i)llo < e + ||w(-,5)||o(is 
provided ||wo|lo ^ ^- 1^ follows that 



If I 



from Lemma 7.7 of Alfaro, Hilhorst, Matano [Ij. The proof is completed. □ 

We have shown in Section 3 that for any initial function uq G C(M; [0, 1]), 
the mild solution uq) of ^ exists for any t > and that u{-,t,uo) G 

C(M; [0, 1]). Now define an operator Q : C(M; [0, 1]) x M+ ^ C(M; [0, 1]) by 

Qt{uo){x) = Q{uo,t){x) = u{x,t,uo), 

where u{x,t,uo) is the solution of Q with initial data uo{x). 

In the following proposition, we will show that Q is a semiflow on the 
space C(M; [0, 1]) with respect to the local uniform topology, in other words, 
with respect to the norm || • ||q. 

Proposition 5.4. Q is a semiflow generated by the solution of equation ([3]) 
in C(IR; [0, 1]) (with respect to local uniform topology) in the following sense: 

(a) Qo{uo) = uo, 
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(b) Qtr+t2 = Qti ° Qt2, 

(c) Q is continuous in {uQ,t). 

Moreover, for any t > 0, {Qt(no) | uq G C(M; [0, 1])} is precompact in C(M; [0, 1]). 

Proof. The properties (a),(b) are obvious. To prove property (c), we only- 
need to show that 

(1) For any given uq, Q{uQ,t) is continuous in t, 

(2) For any given T > 0, the family of maps Qt : C(M; [0, 1]) C{R; [0, 1]), 
< t < T, is uniformly equicontinuous. 

(1) comes from Lemma [3.2[ In Proposition [531 we have shown that property 

(2) holds. This completes the proof that Q is a semiflow. 

The compactness of {Qt(uo) I uq G C(M; [0, 1])} is equivalent to the uni- 
form equicontinuity of the function family. So it can be obtained from 
Proposition 13.31 directly. The proof is complete. 

□ 

In the following lemma, we prove that the semiflow Q is monostable. 

Lemma 5.5. For continuous periodic initial data 1 > uq > 0, uq ^ 0, we 
have u{x,t,uo) — )• 1, as t —)• oo. 

Proof. We consider interval [0,L]. We have 

ut = Uxx + b{x)u{l - u) 
with initial data n(x,0) = uq{x). Let v be the solution of 

Vt = Vxx 

with initial data v{x, 0) = uq{x). From the classical theory of the heat equa- 
tion, v{x, t) > for any t > 0, x G [0, L]. 

From the weak comparison principle, we have 1 > u{x,t,UQ) > 0. More- 
over 6 > 0. Hence 

(U - V)t >{U- V)xx, 

which means u{x, t) > v{x, t) > 0. For a given ti > 0, we may find a constant 
a > such that u{x,ti) > a. This can be done because u{x,ti) is periodic 
with respect to x. 

Let v{x,t) be the solution of 

Vt = Vxx + h{x)v{l - v) 

with initial data i)(x,0) = a. Obviously 

u{x, t + ti, no) > v{x, t). 
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Since a > is a sub-solution of equation vt = Vxx + b{x)v{l — v), from the 
comparison principle, we can get that v{x,t) is increasing in t. 
Set 

P'^{x) = lim v{x,tn)- 

Then 

P+ +6(:r)P+(l-P+) = 

in the weak sense. From P^. < and P+(x) is periodic and P^{x) > a > 0, 
we have = 1. Since u{x,t + ti,uo) > v{x, t), we have u{x, t, uq) — )• 1, t — )• 
oo. □ 

Summarizing, for any t > 0, Qt has the following properties: 

(i) Qt is order preserving in the sense that if uo,vo E C(M; [0,1]) and 
uo{x) < vo{x) on M, Qt{uo){x) < Qt{vo){x) on M. 

(ii) QtiTiiuQ)) = TL(Qf (no)), where Tf, is a shift operator with Ti(n)(3;) = 
u{x — L) . 

(iii) Qt{0) = and Qt{l) = 1. For any uo G C(M; [0, 1]) with uo{x+L) = uq 
on R and u ^ 0, Qt{uo) — )■ 1 in the space C(M; [0, 1]) with respect to 
the local uniform topology. 

(iv) Given T > 0, the family of maps Qt : C(M; [0, 1] ^ C(M; [0, 1])), < 
t < T, is uniformly equicontinuous with respect to the local uniform 
topology. 

(v) Qt{C{M.; [0, 1])) is precompact in C(M; [0, 1]) with respect to the local 
uniform topology. 

Thanks to these properties, the theorem holds: 

Theorem 5.6. For any b G A(a), the spreading speed in the positive and 
negative directions c**{b) and c**{b) exist and they coincide with the minimal 
travelling wave speed in the positive and negative directions c*{b) and c*{b) 
respectively. 

Proof. This theorem can be obtained from the results of Weinberger [H] or 
a more abstract results of Liang, Zhao [8]. In fact, for any t > 0, Qt satisfies 
(i)-(v) which are Hypotheses 2.1 in [l3], and from Theorem 2.6 of p!3] we 
can get the existence of c*(6), b G A{a). □ 
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5.2 Proof of c*(fe) = c*(6) 

To prove c*{b) = c*(6), we consider the linearized equation 



+ b{x)u. 



We define a linear space X by 



X :={(/> 



1,6 GK, </>i,</.2 G 5C(M)} 



and the subset 




{0 G X 



for any M, ^ > 0. Again we equip X with the local uniform topology too. 
Similar to the proof of Theorem 13.11 and Proposition 15. 3|, we can obtain 

Lemma 5.7. For any G X, the mild solution u{x,t) of the equation 



with initial data u{x, 0) = (j){x) exists for all t > and is unique. The mild 
solution is a weak solution and for any t > 0, u{-,t,(/)) S X. Moreover, the 
mild solutions depend on the initial data continuously with respect the local 
uniformly topology on X^^ for any > 0. 



with initial data n(x,0) = uq{x). Similarly to the proof of Lemma 16.31 and 
Proposition 15.41 we can prove the following lemma: 

Lemma 5.8. For any t > 0, is continuous on X with respect to the local 
uniform topology. 

Next, for any M > 0, let A% = C(M; [M, M]) and vl^^ = {e?^n | u £ A^} 
for G M. Then we have 

Lemma 5.9. $t(j4^) is precompact in X with respect to the local uniform 
topology. 

In the following lemma, we show that is strongly order-preserving. 

Lemma 5.10. For any uq G BC{R) with Uo{x) > 0, Vx G M, n ^ and 
any t > 0, we have $j(uo)(x) > for x G M. 



ut = Uxx + h{x)u 



Define 



^tmx) = 

where n(x, t, uq) is the solution of 



u{x, t, (p), ycj) £ X, 



Ut = Uxx + b{x)u 
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Proof. Let u{x, t, uq) be the solution of 

ut = Uxx + b{x)u, 
with initial data uq, and let v{x,t) be the solution of 

with initial data uq{x). 

Since b(x) > 0, we see that v{x,t) is a sub-solution of ut = Uxx + b{x)u 
with initial data uq{x). Then we get 

u{x, t) > v{x, t). 

From the classical theory of the heat equation, v{x,t) > for i > 0. So we 
have ^tiuo) = u{x,t,uo) > v{x,t) > 0. □ 

Let r = {a{x) G C(M) n L~(M), a{x) = a{x + L)}. We equip T with the 
local uniform topology, which is also equivalent to the L°° topology on F. 
For any ^ > 0, define a linear operator jC^ on F by 

4(a) = e«^$t(e-«^a). 
From the definition and the properties of <l>f, we have 

Lemma 5.11. For any t > and ^ > 0, : F — > F is bounded, compact 
and strongly positive. 

Proposition 5.12. ip is a principal eigenfunction of —L^-^ if and only if il) 
is a principal eigenfunction of C^. /x(A, 6) is a principal eigenvalue ofL^j^ if 
and only i/exp(— /x(A, 6) + A^) is a principal eigenvalue of JC^. Consequently 
the principal eigenvalue of L^ i is unique. 

Proof. Since ^ is a principal eigenfunction of L^ i, 

ip" — 2A^' + b{x)il) = — /x(A, in the weak sense. 

Letcf) = e~^^ip. Then the above formula is equivalent to 

(p" + b{x)(/) = {—fi{X, b) + A^)(/> in the weak sense. 
Change begins It is easy to show that exp((/i(A, 6) + X'^)t)ip is the mild 



solution of Uf = Uxx + bu with uo{x) = iIj{x). Change ends Hence, 
$t(0)=exp((-/i(A,6) + A2)t)</,. 

Finally, we have 

A^(V)=exp((-M(A,6) + A2)t)V. 

Since is a strongly positive compact operator, by the Krein-Rutman 
theory, its principal eigenvalue is unique, hence that of L^ i. This completes 
our proof. □ 
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Remark 5.13. Combining Proposition \5.1^ and Lemma \4. 7[ we know that 
Proposition \2.12\ holds. 

Proposition 5.14. The principal eigenvalue of —L^i is negative, and then 
the principal eigenvalue of is lager than 1 for any t>0. 

Proof. Integrating from to L, we get 

fm ^^ b(x)u(x)dx 

^{0,b) = < 0. 

/o u{x)dx 

Furthermore, the principal eigenvalue of £j is e"'^^'''^)* > 1. This com- 
pletes the proof. □ 

Summarizing, for any t > 0, has the following properties: 

(I) is strongly order-preserving in the sense that for any uq £ i?C(M) 
with uo{x) > 0,Vx G M,n ^ 0, <^t{uo){x) > 0, for xeR. 

(II) <I>t(?L('Uo)) = TLi^t(,uo)), where is a shift operator with Ti{u){x) = 
u{x — L) . 

(Ill) For any t > and ^ > 0, the operator £f : F — )• F is bounded, 
compact and strongly positive. Moreover, the principal eigenvalue of 
is larger than 1. 

In the following lemma, we will show that Q can be dominated by ^ from 
above. 

Lemma 5.15. Let Qt and <I>j be defined as above. Then for any uq with 
uo{x) > 0, we have Qt{uo) < ^t{uo) for any i > 0,x G M. 

Proof. We recall that Qt{uQ){x) is the solution u{x,t,uo) of the equation 

ut = Uxx + b{x)u{l - u) 
with initial data uq and (^t{uQ){x) is the solution v{x,t,UQ) of the equation 

vt = Vxx + Hx)v 

with initial data uq. 

From the weak comparison principle, < u < 1. Let w = u — v. Then 

wt < Wxx + bw. 

From the classical comparison principle of heat equation, we have w < 0. So 
Qt{uo) < ^t{uo) for any t > 0,x G M. □ 
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On the other hand, we want to find some Unear operator dominate Q 
from below. 

Lemma 5.16. For any e with < e < 1, let be an operator such that 

m^o) = u^{x,t,uo) 
where n'^(x,t, mq) is the solution of 

ut = u^xx + (1 - ^)Kx)u^ 

with initial data uq G BC(M). Then <I>^ also has the properties (I-III). 
Moreover, for any given to > 0, <I>^j, (uq) < Qto(^o) provided Uo{x) > and 
max|Mo| 'is small enough. 

Proof. We recah that Q{x,t,uo) is the solution of ut = Uxx + b{x)u{l — u) 
with initial data uo{x) and $f (mq) is the solution of = u^.^ + (1 — e)b{x)u'^ 
with initial data uq. Let w = u — u^. Then 

wt = Wxx + &((1 - u)u - (1 - e)n^). 

Since to and e are given, from the continuous dependency of u{x, t) on the 
initial data, if maxuo is sufficiently small, then u{x,t) < e, for t G [0,to]- 
Then we have 

Wt = Wxx + S((l - u)u - (1 - e)n^) > w^x + ^1 " ^^)^- 
It can be rewritten as 

Wt > Wxx- 

Hence, w = u - > 0, for x £R,t e [0, to]. □ 

Proof of Theorem \2.1b\ To prove this theorem, we use the results of Wein- 
berger ^14j or the more abstract results of Liang, Zhao [8j. In fact. Theorem 
15.61 shows the existence of c*{b) . Moreover, for any t > 0, we can check that 
$t satisfies the hypotheses of Theorem 2.5 in [T3] and for any 1 > e > 0, 
^\ satisfies the hypotheses of Theorem 2.4 in [T^. Moreover, for any A > 0, 
the principal eigenvalue of Uxx — 2Xux + (1 — e)b{x)u under the periodicity 
conditions converges to the eigenvalue of Uxx — 2Xux + b{x)u under the same 
periodicity conditions as e — )• 0. Hence we get 

c*{b) = cl{b), bGA{a). 

Similarly 

c*(b) = cl{b), beA{a). 

To prove c*{b) = c*{b), we first prove that, for smooth b G A(a), c*(6) = 
c*(6). Then by using Proposition 14.71 and the result that c*{b) = c*(6), we 
get c*{b) = c*{b). 

For smooth b G A{a), consider the following two operators: —L\^i,^ = 
—ij;"+2\ip'—bip and —L_x^i,'ijj = —jp"—2X'ip'—btlj. They are adjoint operators, 
so they have the same principal eigenvalue which means /x(A, b) = A, b). 
By Definition it holds that c*(6) = c*(6). □ 
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5.3 Proof of Theorem ICTfl and [2A8] 

Now we are ready to prove Theorem 12.171 and 12.181 
Proof of Theorem \2.17\ By Lemma 14.101 

c*(/i) = sup c*(6). 

6eA(a) 

On the other hand, in the previous subsection, we have proven that c*(6) = 
c*{b) for b G A(a). Since h G A(a), we get 

c*{h) = sup c*{b) = rnax c*{b). 

b£A{a) beA(a) 

□ 

Proof of Theorem \2.18[ Theorem 15.61 imphes that 

c**(b) = c*{b), c**(b) = c*(b). 
By Theorem 12. 16^ c*{b) = c*(b). Consequently, 

c**(b) = c**(b). 

□ 

6 Proof of the lemmas 

In this section, we prove the following technical lemmas on the equicontinuity 
of solutions of Cauchy problem ()14p . These lemmas have been used in 
Sections 3 and 5. 

Lemma 6.1. Let u{x,t,UQ,b) be the solution of ()14p with b replaced by a 
smooth b G A(a) with initial data uq G C(M) nL°°(M). Then for any e, M > 
0, {^i(a;, t, ^)}t>e,||uo||<Af,feeA(a) ^-5 uniformly equicontinuous in x G M. 

Lemma 6.2. Let u{x,t,UQ,b) be the solution of (|14p with b replaced by 
a smooth b G A(a) with initial data uq G C(M) H L°°(R). Then for any 
e, M > 0, {u{x,t,uo,b)}x^M.^\\uo\\<M,b&A{a) uniformly equicontinuous in 
t G [e, oo). 

Prom Lemma 16.11 and Lemma 16. 2^ The following lemma easily follows. 
We omit the proof since it is straightfoward. 

Lemma 6.3. Let u{x,t,uo,b) be the solution of ()14p with b replaced by 
a smooth b G A(a) with initial data uq G C(M) H L°°(R). Then for any 
e, M > 0, t, ■uq, &)}||„(,||<^/ ;,£A(q) uniformly equicontinuous in {x,t) G 

M X [e,oo). 
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6.1 Proof of Lemma 16.11 



Given b and uq, we denote u{x,t,uo,b) by u{x,t) simply. First, we recall 
that u{x, t) is a mild solution. It can be written as 



u{x, t) 



'AlTt 



e « u{y,0)dy 



+ 




\/47r(t - ^ 



--biy)u{y, s){l - u{y, s)) dyds). 



Let t > be given. In order to prove that u{x,t) is equicontinuous with 
respect to x, we need to prove that for any small e > 0, there exists a 5 > 
such that 

\u{xi,t) - u{x2,t)\ < e, 
if |xi — X2I < S. For any xi > X2 £ M, 



u{xi,t) - u{x2,t) 



1 



[e it 



)u{y,0)dy 



+ 



g 4(t-s) g 4(t-s) 

V'4vr(t - s) 



bn{y)u{y, - u{y, s)) dyds. 



Next we prove that for any e > 0, there exists a > 0, such that the 
first part 



and the second part 




g 4(4-3) g 4(t — s) 



JR 



)u{y,0)dy 



bn{y)u{y, s){l - u{y, s)) dyds 



can be bounded by e if |xi — X2I < 6. 
For the first part, 



1 



(e 4t 



1 



-e 4t )u{y,0)dy 

-e 44 u(2/,0)c(y X (xi - X2), 



/Ant Jr 4t 
where function of 1/ satisfying 

(^i-y)^ 2(xc y — y) _(fijr:Zl! 

g 4t — g 4t = — ^^i^ -g 4t 

and x^^y € (x2,xi). We know that the integration 



Ixle ^ dx 
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is bounded. Now we need to prove that 



V ^.2/ 4^ — u{y,0)dy 



At 

is bounded. Consider the following sets: 

= {y G M|i < I sup{|?/ - xi\,\y - X2\} < i + 1}, i = 1,2,3.... 

We have 

1 f 2{x^^y-y) H,y-yf 



At 



'u{y,0)dy 



-1 2{x^^y-y) 



VA^VU VAi 



e 4t u{y,0)dy. 



There exists a constant M{t) > 0, such that for every ilj, 



1 2{x,y-y) _ H,y-y)' 



n-VA^VAi VAi 



u{y,0)\dy 



< 



M{t) [ \y- 



x\\e 



dy 



which means the first part can be bounded by e/2 if \x\ —X2\ is small enough. 
Next step is to prove that the second part 



g i(t-s) _ g 4(t-s) 



Jm. \/A'K{t — s) 



Ky)u{y, s)(l - u{y, s)) dyds\ < e/2. 



if \xi — X2I < 5. 

For notational simplicity, in what follows we write: 

supg{x) := sup g{x). 

x,kL xe[kL,kL+L] 



First, 



4(t-s) 



Jr VAiT{t — s) 



Ky)'<^{y^ - s)) dyds\ 



< M / V supe"^^^ ds =^ M'{t), 

Jo V47r(t - s) ^j/,fcL 

where M > is a constant related to the bound of u and M'{t) is bounded 
for fixed t. Hence 



g 4(t-s) 

Jr VA-Kit - s 



-.b{y)u{y, s)(l - u{y, s)) dyds 
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is uniformly bounded for a; G M. Then we can choose at*<t such that 



e 4(t-s) 

/ /A M y)u{y, - u{y, s)) dyds\ < e/4 

t* JR \'4:TT{t - S) 

Ky)^{y^ - ^{y^ «)) dyds\ < e/4, 



and 



It* JR \/47r(i — s) 

where t* is independent of the choice of xi and X2- 
Consider the integral from to i*. 



/■** f e 4(t-3) _ e 4(t-s) 

/ / ,^ b{y)u{y,s){l -u{y,s))dyds 

Jo Jm. \/4:Tr[t — s) 



■b{y)u{y,s){l-u{y,s))dyds{xi-X2) 



10 J-R ^47r(t — s) 

= r / 

Jo J^irit - s)J4:{t - s) Jr 



, , , ^* b(y)u(y, s)(l—u(y, s)) dvds 

V47r(t-s)V4(t^yR V47r(t - s) " ^ 

•x{xi - X2), 

where x^^y is a function of y and takes values between X2 and xi. 
Since t* < t, we can get 

' -Ky)uiy,s)il-u{y,s))dyds\ 



r -\ f 

Jo A/47r(t^^v/4(t-s) J^ 



y/A7r{t- s)^/4{t- s) Jr ^/4Tr{t - s) 
is bounded. 

Then if we choose S small enough, when |xi — X2\ < S, we have 



Jo VmK{t - s) Jr 



s) 

X (xi -X2)\ < e/2. 



-b{y)u{y, s)(l - s)) dyds 



It is not difficult to see that 6 is independent of xi and X2. 

Combining the above estimates, we see that there exists (5 > such that 
if — X2\ < S, then \vn(xi,t) — Vn{x2,t)\ < e. It is not difficult to see that 
the above estimate is independent of the choice of b and uq. This completes 
the proof of the lemma. 
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6.2 Proof of Lemma 16.21 

As in lemma EH we still use 



u{x,t) = —^= e 4t u{y,0)dy 
\/47rt Jr 

—?===Hy)u{y, s){l - u{y, s)) dyds. 

Here we consider an interval [Ti,T2] C M, where T2 > Ti > 0. For 
ti<t2e [TuT2], 

u{x,tl) - U{x,t2) 

1 f (^-i/)^ 1 f (^-y)^ 

e **! u{y, 0)dy -= / e **2 0)(iy 



e 4(ti-s) 

+ / / /. / =b{y)u{y,s){l-u{y,s))dyds 
\/47r(ti - s) 




/o 

ng"4('t2°'s) 

We divide the right-hand side into two parts as in Lemma |6.1[ 
1 f (^-i/)^ 1 /■ (^-i/)^ 



and 



{I):=^==J^e uiy,0)dy - ^== I e ^'2u{y,Q)dy 



ne 4(ti-'') 
/. M y)u{y, - ^/(y, s)) dyds 

t2 f g 4(t2-s) 

b{y)u{y, s){l - u{y, s)) dyds. 



'0 



\/4vr(i2 - s) 



First, 

(/)= / — , — -(ti-i2My,o)dy, 

where i^^^y is between t2 and ti. 

As discussed in Lemma 16.11 we may get that (/) can be very small if 
1^1—^2! is small uniformly in G R. 
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The next step is to prove that for any e > 0, there exists a5 > such that 
the absolute value of the term (//) is less than or equal to e if \ti — ^2! < 



Jo JR 



Jr Y^47r(ti - s) 



biy)uis, y)(l - u{y, s)) dyds 



rt2 




e ■4(*2-'<) 




V^47r(t2 - s 



Myhiy, s){l - u{y, s)) dyds 



_(x~yf_ 



-My)u{y^ ■s)!! - u{y, s)) dyds 



Y^47r(f2 - s) 



ft2 



g 4(t2-s) 

Iti Jr •y/47r(t2 - s) 
It is not difficult to show that 




biy)uiy, s){l - u{y, s)) dyds 



Ky)u{y, s){l - u{y, s)) dyds. 



rt'z 




g 4(t2-s) 



--Ky)u{y,s){l-u{y,s))dyds < D^/\ti - t2\, 



^/Mt2 - s) 

where D is a constant. Therefore it suffices to consider 




g 4(ti-s) 

V47r(ti - s 

tl 



--b{y)u{y, s){l - u{y, s)) dyds 




_{x-yf_ 
g 4(t2-s) 



-Myhiy, s){l - u{y, s)) dyds. 



'0 Jr v^47r(t2 - s) 
Once again, as in lemma [6Tl for e > 0, there exists a t* < ti, such that 



g 4(ti-s) 

t* Jr y^47r(ti - s) 



b{y)u{y, - u{y, s)) dyds\ < e/4 



and 



g 4(t2-s) 



2-S 



--b{y)uiy, - u{y, s)) dyds\ < e/4. 
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The left-hand side is an integral from to t* which is 

e 4(ti-s) g 4(t2-s) 

( /, ,^ ^ - /, . )Ky)u{y, - u{y, s)) dyds 




ft* 

10 




X {ti - t2)b{y)u{y, s)(l - u{y, s)) dyds 

* -{ti-t2)b{y)u{y,s){l -u{y,s))dyds. 



Next we prove that both 



4(ts,j, s) y/^^m,y s) b{y)u{y,s){l -u{y,s))dyds 



Jr ^^ik,y - 





'0 



and 



b{y)u{y,s){l-u{y,s))dyds 

are uniformly bounded in a: G M. Since they are similar, we prove the former. 
We have 

(x-y)'^ 

<M| r ^ ^, ^ , , / {x-yfe~^^b{y)dyds\ 
As we have proven above, 

/ {x-yfe~^*i^^b{y)dy 
Jr 

is uniformly bounded for any a; € M. Since t* <ti, this integral is bounded 
by a constant M{t*) > 0. Consequently there exists a di > 0, such that if 
t2 — ti < Si, then 

4(t,„-s)VM*e..-^) (ti-t2)6(yMy,5)(l-u(y,.))dyd5| < e 
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We can also prove that there exists 62 > such that if t2 — h < 62, then 



Here ti,t2 G [Ti,T2] and the equicontinuity is independent of the choice of 
X, uq and b. The proof of the lemma is complete. 
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